Island phases and charge order in two-dimensional manganites 
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The ferromagnetic Kondo lattice model with an antiferromagnetic interaction between localized 
spins is a minimal description of the competing kinetic (t) and magnetic (K) energy terms which 
generate the rich physics of manganite systems. Motivated by the discovery in one dimension of 
homogeneous "island phases", we consider the possibility of analogous phases in higher dimensions. 
We characterize the phases present at commensurate fillings, and consider in detail the effects of 
phase separation in all filling and parameter regimes. We deduce that island and flux phases are 
stable for intermediate values of K/t at the commensurate fillings n — 1/4, 1/3, 3/8, and 1/2. 
We discuss the connection of these results to the charge and magnetic ordering observed in a wide 
variety of manganite compounds. 



PACS numbers; 75.10.-b, 75.30.Vn, 75.40.Mg 
I. INTRODUCTION 

Transition-metal manganite compounds have long 
been known to display a broad spectrum of physical prop- 
erties as a function of temperature, filling and counte- 
rion composition. While the most remarkable of these 
is the colossal magnetoresistanceld observed in the ferro- 
magnetic (FM) phase, the phase diagrams of both cubic 
perovskite and layered manganite materials exhibit a rich 
variety of metallic, insulating, magnetically ordered and, 
apparently, inhomogeneous or phase-separated regions. 

The ferromagnetic Kondo lattice model (FKLM) has 
been used extensively as a minimal model to reproduce 
the physics responsible for this situation. We will study 
a version of the model which includes a Heisenberg in- 
teraction between the localized spins. In essence, this 
encapsulates the competition between the ferromagnetLe 
polarizing effect of the double-exchange hopping terniEl 
(t) for mobile carriers in the Cg orbitals of Mn'^"'", and 
the antiferromagnetic (AF) interaction (K) between the 
localized spins composed of electrons in the t2g orbitals. 
Treatments of the model with both classical local spins, 
and with fully quantum, 5 = 1/2 local spins, both re- 
turn some of the features observed among the selection 
of manganite phase diagrams. A large number of authors 
has worked on many forms of the FKLM, and we will 
present in the following sections only a small selection of 
references relevant to the current approach. 

Following the discoverjrl in one-dimensional simula- 
tions of novel "island phases" near commensurate val- 
ues of electron filling in the FKLM with strong Hund 
coupling between localized and conduction electrons, 
we wish here to consider the possibility of higher- 
dimensional generalizations of these phases. By an is- 
land phase is meant a spin configuration composed of 



small, regularly arranged, FM islands (clusters of 2-4 
sites in Ref. ^), with AF local spin orientations between 
islands (Fig. 1). These phases are homogeneous, and near 
the commensurate fillings maximize kinetic energy within 
each island at minimal cost to the magnetic energy, which 
is favored at the island boundaries. Focussing primarily 
on the problem in two dimensions (2d), we wish to estab- 
lish the possibility that such islands, which may be small 
in one or both directions, remain the most stable phase 
for certain fillings and parameter ratios K/t. 
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FIG. 1. Schematic representations of the island phases 
(7r/3,7r) (a) and (7r/2,7r/2) (b). 

A particular motivation for our study is the recent 
observation of charge-ordering phenomena, and more 
general inhomogeneous charge and spin configurations, 
in a variety of manganite systems. These appear in 
both layered and cubic materials, and at both commen- 
surate and incommensurate values of the electron fill- 
ing set by the counterion doping. Some of the ear- 
liest observations of charge orderinga were made in 
Lai„a;Sra;Mn03, and were followed bji-measurements sug- 
gesting ]pplarons,El phase separation,^ and paired stripe 
features.!] Charge order coupled to a structural phase 
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transition has been observed in Bai-jjCjiajMnOa at in- 
commensurate values of the filhng xJaa Among hole- 
doped manganites, charge ordering arose at incommen- 
surate filling in Ndi-j^Sr^MnOs, and in a stppj&ilike con- 
figuration at half-filling in Pro.5Sro.5MnO3.By For the 
latter system, the stripe feaijures could be made to "melt" 
in an applied magnetic field. t3 Of most interest in the cur- 
rent context, ordering phenomena have also appeared in 
2d or layered mangapije systems. In Sr2_2;Laa:Mn03 at 
low doping, Bao et aZjlj reported charge order, phase sep- 
aration, and triplet bip(j>larons. For the same system at 
X — 0.5, Moritomo et a/.t2l related charge ordering to laU 
tice effects by substitution for La, and Murakami et alM 
made direct measurements of charge and orbital order 
for the commensurate La member. Finally, we mention 
also the observationEj of charge order in the layered 327 
compound LaSr2Mn207. 

The manuscript is organized as follows. In Sec. II we 
present the model in the form we wish to consider, and 
outline the methods by which it is analyzed. In Sec. Ill 
we discuss the available means to characterize the phases 
which appear, and illustrate these with examples. Sec. 
IV contains a detailed discussion of the issue of phase sep- 
aration, and a global phase diagram for the augmented 
FKLM which delimits the regimes of interest for island 
phases. We return in Sec. V to the robust fiux and island 
phases, discuss their properties and their charge order, 
and consider their relevance to the above experiments. 
Sec. VI gives a summary and conclusions. 

II. MODEL AND METHODS 

We consider the FKLM in the form 

H = - ^ tij (cl^Cia + H.c.^ - Si. Si 

{ij)a i 

+if5]S,.S,. (1) 

(ij) 

Here cj^ is the operator creating an electron of spin a 
in the sole Cg orbital; = X]a/3 '4a'^a(3Cii3 gives the spin 
of this "conduction" electron, and its mobility depends 
on the orientation of the localized t2g spins according to 
the double-exchange mechanism.0 The second term is the 
Hund coupling, Jh > 0, which favors a FM orientation 
of spins on the same site. Following Refs. |l6|,0 we will 
be concerned with the limit of large Jh', while in real 
systems Jh is of the same order as the bandwidth, this 
simplifying approximation has been found to give reason- 
able results. The limit corresponds to a situation where 
the conduction electron is bound to follow the spin tex- 
ture of the localized system, while anti-aligned electrons 
occupy a band with energy higher by Jh- The projecting 
effect of the large Hund coupling allows one to neglect di- 
rect Coulomb interactions of the eg electrons. The final 
term, with K > 0, expresses the AF interactions between 



the local t2g spins, whose competition with the FM spin 
alignment required to maximize Cg electron kinetic en- 
ergy (H) generates the intrinsic physics of interest in the 
context of manganite materials. The situation is repre- 
sented schematically in Fig. 2. 
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FIG. 2. Schematic representation of the Hamiltonian (1) 
for two sites. 

We will analyze the model primarily by a classical 
Monte Carlo (MC) procedure for the localized spins, in 
conjunction with-esact diagonalization of the conduction 
electron system.tMj The localized spins are thus taken 
to be classical, an approximation to the true situation of 
5 = 3/2 which is found not to invalidate the connection 
to real systems. The conduction electrons are taken to 
occupy a single eg orbital, or band, and from the condi- 
tion on Jh only one spin projection need be considered. 
This part of the process is the solution of the single- 
electron problem with hopping set consistently by the 
localized spip configuration. In the limits of large S and 
Jh, this iM 

t ^ t f cos - cos ^ + e-('^--'^.) sin ^ sin ^-A , (2) 
^ V 2 2 2 2/^^ 

where 9i and (pi are the polar angles of spin S.^. The 
resulting energy levels are then filled by the available 
number of electrons in the canonical ensemble. 

The MC simulation proceeds from the FKLM partition 
function with classical spins, 

N-kN 2n 

Z=T\ de^sinO^ Tr[exp(-/3F)], (3) 

,^1 Jo Jo 

where N is the system dimension. Positivity of the inte- 
grand assures that the sign problem is absent. Updates 
of the spin configuration {0i, (pi} are accepted or rejected 
according to the Glauber algorithm. Because there are 
no cases where non-coplanar spin configurations appear, 
and because of the large degeneracy of coplanar phases, 
the simulations could be accelerated by fixing 9i = tt/2, 
and varying only the angles {(pi}. The number of MC 
steps per site for iV = 8 is taken as 2000 to equilib- 
rium and 3000 for measurement, while for = 12 the 
corresponding numbers are 500 and 1000. Systems of 
size up to 12x12 are accessible by this method, and thus 
we supplement the MC results by a variety of classical 
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analytical considerations, which afford considerable in- 
sight and allow a detailed assessment of finite-size effects. 
The simulations may be pursued down to temperatures 
of T = 0.005t, which unless otherwise stated will be the 
relevant value for the MC results displayed. This temper- 
ature is sufficiently low that comparison with the zero- 
temperature, analytical calculations is meaningful, and 
in most cases quantitatively so. The method is the same 
as that used hsL Dagotto, Yunoki and coworkers in a series 
of papers JIJlj We will reproduce some of the same re- 
sults, and comment on the similarities and differences in 
the context of our island phase analyses in what follows. 

III. PHASE CHARACTERIZATION 

In this section we will present some results for typical 
phases which emerge from MC simulations performed at 
the commensurate fillings n — 1/2, 1/3, and 1/4, and 
for the full range of values of K/t. The results of the 
simulations for the localized spin system may be charac- 
terized by three separate but related quantities: the spin 
structure factor 

5(k) =^S,.S,e*-('-'-'-^), (4) 

a histogram of the distribution of angles between all 
nearest- neighbor spin pairs, which we choose to present 
as a function of cos Qij , and a simple "snapshot" of the 
spin configurations at a representative step late in the 
MC process. Note for the histogram that 6^^ is the full 
angle between spins given by cos = (Sj.Sj)/^^ for the 
classical case, and is not to be confused with the on-site 
azimuthal angle 6i in Eq. (^. Finally, one may compute 
in addition the charge distribution function 

n(k)=^n,e'''--, (5) 

i 

and (by analogy with Eq. (^)) the charge-charge corre- 
lation function N(k), which we will use in Sec. V when 
considering charge order. 

As a guide to understand the variety of possibilities 
which is contained in these quantities, we first calculate 
the classical, ground-state energies of a multiplicity of 
possible spin configurations. This may be carried out for 
an infinite 2d system by straightforward extension from 
the arguments presented for the Id case in Ref. |. For 
each spin configuration, the magnetic energy per spin is 
a simple function of the average of the angles across each 
bond, which varies from 2K for the FM case to — 2/\ 
for the AF. The kinetic energy at this level is a readily 
calculable function of the spin configuration which varies 
from in the AF case, where all kinetic processes are ex- 
cluded, to the average energy of the 2d nearest-neighbor 
band = —2t (cos + cos ky) , for the relevant band 
filling, in the FM case where it is maximally negative. 



The results of this exercise are illustrated in Fig. 3 for 
n = 1/2, in Fig. 4 for n = 1/3, and in Fig. 5 for n = 1/4. 

All of the phases denoted by {kir/m, In/m) have neigh- 
boring spins only either parallel or antiparallel, in both 
directions. The rational fractions k/m,l/m may be un- 
derstood as indicating that the spin direction turns over 
k or / times in 2m lattice constants. Fig. 1 shows two 
small-rn possibilities, the {tt/3,tt) (a) and (7r/2,7r/2) (b) 
phases. As a more complex example, the phase (37r/4, tt), 
which appears over a wide range of K/t at filling n = 1/4 
(Figs. 5,8), would be composed of chains with repeat unit 
TTitiiti in the x-direction, and AF alignment in the y- 
direction. In addition to these phases, which include the 
FM (0, 0) and-AF (tt, tt) endpoints, we include also the 
"flux phase" £Scj which will be discussed in more detail 
below, and a "double spiral" (DS) phase, by which is 
meant a single phase where the nearest-neighbor spins 
rotate by the same angle < < tt in both x and y 
directions. In this last case, the optimal angle Q is ob- 
tained by minimizing a function of K/t, and the double 
spiral may be expected to be more favorable than any 
variety of single-spiral phases combined with other forms 
of modulation in the transverse direction. Although we 
have considered many possible phases of the above types, 
in Figs. 3-5 we include for clarity only those which are 
the ground state for some range of K/t. 




K/t 

FIG. 3. Energies of selected spin configurations for n = 1/2 
at all values of K/t. Note the clear succession of the ground 
state with increasing K/t from flux phase to (7r/2,7r) to 
{n/2,n)+{n/,-K). 

The calculation of all of these phase energies is 
straightforward. In brief, calculation of the only 2d band 
at (0, 0) proceeds as above, with the filling determining 
the chemical potential up to which the filled band is in- 
tegrated. For the Id structures (0, ^7r/m), one may con- 
sider the band — ~2t cos k in the continuous direction, 
split appropriately into 2, 3, or 4 (the maximum included 
here) by an equal interchain hopping t. Integration over 
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the filled parts of these bands up to the chemical po- 
tential yields the average kinetic energy. For the "Od" 
structures (fc7r/m, Z7r/m), the kinetic energy is a simple 
m^/fcZ-site diagonalization problem to obtain the discrete 
levels. These phases are particularly favorable when the 
filling exactly matches a large gap in the few-level spec- 
trum, e.g. (7r/3,7r) for n = 1/3 (Fig. 1(a)), or (7r/2,7r/2) 
for n — 1/4 (Fig. 1(b)). The calculation of the kinetic 
part for the double-spiral phase follows the 2d case above, 
with reduction of the bandwidth by a factor of cos 0/2, 
while the magnetic part varies as cosG. We do not find 
that canted states are favored in these considerations. Fi- 
nally, two special configurations which require separate 
consideration are the (7r/2, 7r)-|-(7r, tt) phase, to which we 
return in Fig. 14, and flux phases. 
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FIG. 4. Energies of selected spin configurations for n — 1/3 
at all values of K/t. Note the competition of several phases 
around K/t = 0.1. 

Flux phasesH are an important feature of the model 
in any dimension higher than 1. From Eq. ^ it is clear 
that the hopping term also contains a phase factor, and 
that for certain spin textures this phase may differ de- 
pending on the path through the lattice taken between 
two points. The simplest flux phase is that appearing 
at half-filling over a broad range of intermediate K/t, as 
discussed in Ref. |2^, and shown in the snapshot in Fig. 
13(c) below. The term "fiux phase" is used here to refer 
to any spin configuration with this non-trivial topologi- 
cal property, which can be quantified by a non-zero spin 
current .E3 In principle, a variety of fiux phases may exist, 
but we have not yet been able to find any others which 
are ground states at any filling. At the analytical level, 
the semimetallic density of statesc3 of the dispersion 



± \ I cos k1 + cos fc^ 



(6) 



of the simplest fiux phase, which is zero precisely at half- 
filling, accounts for its particularly low energy at rt = 1/2. 



We will characterize this phase in detail in Sec. V. 

While these classical, zero-temperature pictures turn 
out to be rather valuable, and also not quantitatively 
unreasonable, for understanding the 2d pictures to fol- 
low, they arc limited by the imagination of the authors 
as further possibilities may not be excluded. We have ob- 
tained many of the phases proposed in Figs. 3-5 in MC 
simulations, and the following Figs. 6-8 illustrate some 
representative results. 
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FIG. 5. Energies of selected spin configurations for n = 1/4 
at all values of K/t. Note the dominance of the phase 
(7r/2,7r/2) at intermediate K/t. 

In Fig. 6 is shown S'(k), histogram and snapshot infor- 
mation for a phase at filling n = 1/2 and for the ratio 
K/t = 0.22. We see a single peak in S{k) (Fig. 6(a)) only 
at (7r/2,7r), indicating an island phase of FM pairs (the 
"islands") arranged in an AF pattern. The histogram 
(Fig. 6(b)) shows essentially only angles of and tt, rul- 
ing out a possible interpretation as a 7r/2 spiral in one 
direction; the ratios of angles to angles tt is approxi- 
mately 1:3 as expected. Finally, the instantaneous spin 
configuration in Fig. 6(c) illustrates that the simulation 
has in fact converged quite well to the expected phase. 
Comparison with Fig. 3 indicates that for the 2d case, the 
value of K/t for a robust (7r/2, tt) phase is that expected 
from the infinite system at T = 0. 

Fig. 7 illustrates the same quantities for filling n = 1/3 
and K/t = 0.25. For this relatively large parameter ra- 
tio, the dominant (27r/3, tt) phase in 5'(k) (Fig. 7(a)) co 
sists of AF chains with spin configuration TTiTTJ-TTi • 
This is one of the primary types of island phase which 
we will mention again in Sees. IV and V. All three parts 
of Fig. 7 show in addition that this phase is not pure 
in the small-system MC simulation, with spin misalign- 
ments across the cluster manifest in a broadening of the 
histogram peaks and residual components in >S'(k). How- 
ever, Fig. 7(b) does demonstrate the absence of interme- 
diate angles which might be expected from any kind of 
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spiral phase. 
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FIG. 6. MC phase for n = 1/2 at K/t = 0.22, calculated 
for an 8x8 system, (a) Structure factor, (b) Angle histogram, 
(c) Configuration snapshot. 

Fig. 8 characterizes the phase arising for n = 1/4 at 
K/t = 0.20. From Fig. 5 we expect the phase (37r/4, tt) as 
ground state, and indeed this is the dominant component 
in S'(k) (Fig. 8(a)). The rather stronger admixture of 
other components arises because the chosen value of K/t 
is close to a phase crossover, and so other possible 8x8 
phases are not entirely absent. These are not reflected in 
the histogram (Fig. 8(b)) because all the pure phases 
present have angles of only or tt, but the snapshot 
(Fig. 8(c)) does show a small amount of misalignment 
between the predominantly AF-oriented spins. We note 
that the expected pure configuration (see below Fig. 3) 
remains rather hard to observe in Fig. 8(c), and ascribe 
this to the mixing problem, and to the effects of fluc- 
tuations on the small-cluster MC calculation. This ex- 



ample illustrates both the need for careful consideration 
of finite-size effects, and the fact that for all commensu- 
rate fillings there exist regions of K/t (close to the line- 
crossings in Figs. 3-5) where the MC results show strong 
mixtures of different phases. We note in passing that for 
all fillings we find pure FM phases at small but finite K/t 
ratios, in accord with zero-temperature, infinite-system 
expectations based on Figs. 3-5. These straightforward 
cases are not shown here. At large values of K/t, small- 
cluster calculations are unable to access the double spi- 
ral phase, and show instead the AF. We defer a more 
detailed characterization of the most interesting phases 
in these figures, namely the flux phase at n = 1/2, the 
(tt/S, tt) phase at n ~ 1/3, and the (7r/2,7r/2) phase at 
n = 1/4, until Sec. V, after addressing the question of 
phase separation. 

The results of Fig. 6-8 were obtained for small sys- 
tems, where finite-size effects are of paramount impor- 
tance. For fillings n = 1/2 and n = 1/4 we may com- 
pare 4x4 with 8x8 MC results, and for n = 1/3 6x6 
with 12x12. These comparisons give already a good 
indication of where, for example, certain of the many 
possible phases are anomalously favored by the location 
of the chemical potential relative to a gap between sets 
of degenerate states. Even more valuable information is 
provided by comparison with the infinite-system results: 
these may be augmented by performing the same calcula- 
tion, placing spins in a fixed configuration and deducing 
the magnetic and kinetic energies, for the system sizes 
4x4 to 12x12 of the simulations (and further for 16x16). 
An effective calibration of the MC results is then possi- 
ble, by which is meant a renormalization to account for 
effects arising only from system size, which is particularly 
important in discussing phase transitions (Sec. V). 

On these finite systems we are unable to observe phase 
transitions, which are replaced by crossovers occurring in 
a finite range of K/t. As we will show in Sec. V, however, 
a certain amount of care is required in interpreting two- 
peak features in 5'(k), because some robust, single phases 
arising at particular values of filling and K/t do indeed 
have more than one characteristic wave vector in small 
systems. Another feature requiring particular attention 
is the possibility of large-unit-cell phases, which cannot 
be accessed in the MC simulations. An example already 
mentioned is the double spiral, which is expected from 
Figs. 3-5 to be the most favorable phase on approaching 
the AF limit, into which this phase in fact passes contin- 
uously. However, at intermediate to large values of K/t 
we must also consider a competing, large-unit-cell (large- 
m) phase of the type ((m — k)'K/m, tt), k <^ m, with only 
and TT angles between the spins, in which the kinetic 
energy gain comes from spins shared between rare FM 
pairs in an otherwise AF structure. These phases are 
compared in the next section. 
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FIG. 7. MC phase for n = 1/3 at if/t = 0.25, calcu- 
lated for a 12x12 system, (a) Structure factor, (b) Angle 
histogram, (c) Configuration snapshot. 

To conclude this section, wc find that island-Hke phases 
are quite ubiquitous at all intermediate values of K/t 
(Figs. 3-5). The FM islands may be restricted in one di- 
rection, giving rise to stripe-hke features, or in both to 
give true islands, depending on the filling. These states 
are also accompanied by flux phases, of non-trivial spin 
texture, in certain parameter regimes. These novel, ho- 
mogeneous phases arise only as a result of the competi- 
tion between the first and last terms in Eq. (0), without 
recourse to additional physics (a discussion of which is 
deferred to a later section). However, we have worked in 
a canonical ensemble and considered only the energy of 
the emerging phases at zero or the lowest temperatures. 
We now turn to the question of phase separation within 
the model. 
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FIG. 8. MC phase for n = 1/4 at K/t = 0.20, calculated 
for an 8x8 system, (a) Structure factor, (b) Angle histogram, 
(c) Configuration snapshot. 



IV. PHASE SEPARATION 

In the previous section we have considered a canonical 
ensemble, meaning fixed particle number, and deduced 
the ground states on the basis of minimal internal en- 
ergy (or free energy at very low temperature). To ensure 
the global stability of these phases we must consider the 
possibility of their separation into regions of distinct and 
different filling. This propensity has been shown in the 
same model applied in ld,t3 by working in a grand canon- 
ical ensemble and observing discontinuities in filling on 
varying the chemical potential. Here we choose to charac- 
terize phase separation from the energy in the canonical 
ensemble, by observing the curvature of this quantity as 
a function of filling. In Figs. 9-11 are shown the energy 
for fillings n between and 1/2, at low, intermediate 
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and higher values of K/t. We note that the energy is a 
symmetrical function for l/2<n<lby electron-hole 
transformation, and do not comment further on this re- 
gion. In these figures are included data from 12x12 and 
16x16 systems, and infinite-system values for the flux 
and double spiral phases. 

In Fig. 9 we see a convex (up) region at low filling, the 
implication of which is a preference for phase separation 
into two regions, one of zero hole content and the other 
whose filling n is given by a Maxwell construction using 
the tangent to the concave part of the curve. The empty 
region would have AF spin configuration, while for this 
low value of K/t the partially filled region would be FM. 
This result confirms that phase separation is an impor- 
tant property of the model, and agrees qualitatively with 
Ref. |l9|. 
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FIG. 9. Energy as a function of filling at fixed K/t = 0.04 
for a variety of phases. The tangent to the curve indicates 
the regime of phase separation by the Maxwell construction. 

For intermediate K/t (Fig. 10) the situation is more 
complex. The convex regime extends over a much 
broader range of filling, but the "curve" is much less 
smooth, as a result of the particularly favorable island 
phases which can be established at the commensurate 
fillings. In fact. Maxwell constructions applied to Fig. 10 
yield for this value of K/t a separation only into phases 
n — Q and 1/4, or into n = 1/4 and n close to 1/2. In- 
spection of Figs. 3-5 shows that for K/t = 0.12 the phases 
at these two fillings are particularly robust, whereas at 
n = 1/3 a crossing between two phases occurs; by con- 
trast, an n = 1/3 phase would be expected as an endpoint 
of such separation for K/t = 0.15, and indeed emerges 
(Fig. 12, below). This result highlights the dominant 
role of the commensurately filled phases, and suggests 
both "high-contrast" and "low-contrast" phase separa- 
tion. By this is meant in the former case the abrupt 
split into zero- and part-filled regions, and in the latter a 
finer phase separation for certain K/t where incommen- 



surate fillings 1/4 < n < 1/2 may undergo separation 
into regions with closely neighboring, more commensu- 
rate fillings. These statements are made systematic in 
the summary phase diagram presented as Fig. 12. 
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FIG. 10. Energy as a function of filling at fixed K/t — 0.12 
for a variety of phases. The solid lines are Maxwell construc- 
tions. 

At large K/t (Fig. 11) the picture changes again. Here 
the finite-system points for commensurate phases show 
the intriguing feature of lying on a straight line connect- 
ing zero- and 1/2-filling. These are the ((m — A:)7r/m,7r) 
phases introduced above, for those values of m small 
enough for the unit cell to fit within the system studied. 
Simple consideration of fixed spin configurations suggests 
that, in principle, phases of arbitrarily large unit-cell size 
are possible, and their energies will fall on the same line. 
From above, the nature of these phases is an AF config- 
uration of spin chains with k up-spin and k down-spin 
pairs contained in an otherwise AF system with unit-cell 
size 2m. In a fully classical system there would be no 
phase separation with filling in the thermodynamic limit 
at large K/t, but instead a continuous evolution of the 
unit-cell dimension to accommodate the added charges. 
In fact the values of k and m are fixed rather simply by 
the filling n, because the phases of this type appearing 
as the ground state are ((1 — n)^, tt), and their energy is 
given from the number of FM pairs and AF bonds as 



E = -2K + n{K - t) 



(7) 



per site. For the commensurate fillings n = 1/ni = 1/2, 
1/3, and 1/4, we recover the island phases of Figs. 3-5. 
These phases appear to have been overlooked in Ref. |l^, 
although the authors were little concerned with the high- 
K regime. 

Our conclusions are summarized in the global phase 
diagram of Fig. 12. The properties of the minimal form 
of the FKLM (Eq. (0)) faU broadly into four regions, 
determined largely by the ratio K/t of the super- and 
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double-exchange energy scales. For the lowest values of 
K/t^ the system separates into AF and FM phases. For 
small to intermediate ratios, 0.08 < K/t < 0.2, there is 
large-scale phase separation into only the island phases 
appearing at the commensurate fillings n = 1/4, 1/3, 
3/8, and 1/2. An exception here is the flux phase, which 
occupies a finite doping region close to n = 1/2. Wc note 
in passing that within our classical formulation, only the 
FM and flux phases offer the possibility of hopping of 
conduction electrons throughout the system; only these 
phases would have metallic properties, and all others will 
be insulating. 
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FIG. 11. Energy as a function of filling at fixed K/t = 0.24 
for a variety of phases. The solid line is a guide to the eye. 

For intermediate ratios 0.2 < K/t < 0.28 we find the 
large-unit-cell phases discussed above. The hierarchy of 
possible states exists across the full doping range only 
when no competing phase falls below the straight-line 
energy function (Fig. 11, Eq. (Q)) for any filling, and it is 
this condition which sets the limits in K/t of the shaded 
region in Fig. 12. We have marked (vertical dashed lines) 
the small-m phases which are compatible with the fi- 
nite clusters considered, but stress again that from the 
present calculations we expect to find all phases of the 
form ((1 — n)TT,Tr) for the infinite system. All states in 
the shaded region are a form of 2-site FM island phase, 
which would show charge-ordering peaks in Af(k) (Sec. 
V), while the small-m members at the commensurate fill- 
ings provide examples which may be studied on small 
clusters (Figs. 6-8). At intermediate to large values of 
K/t, the large- unit-cell phases are replaced by a wide 
region of "high-contrast" phase separation due to the 
extraordinary stability of the (7r/2, 7r)-|-(7r, tt) phase at 
n = 1/2. We have found only this phase, which is con- 
sidered in more detail in the following section, and the 
AF phase with zero filling, to be stable in this regime 
of K/t,hut stress that we cannot fully exclude the possi- 
bility of similar ((m — /c)7r/m, 7r)-|-(7r, tt) phases at other 



commensurate fillings. A search for these is limited by 
the available cluster size, and remains a topic for future 
investigation. Finally, at large values of K/t we recover 
the conventional, spiral-ordered DS phase, which passes 
smoothly to an AF phase. 
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FIG. 12. Phase diagram of augmented FKLM for the full 
range of filling n and ratio K/t. PS denotes phase separation, 
the thick, vertical lines the island phases, and the shaded 
region the regime of large-unit-cell phases. 



V. ISLAND PHASES 

With the results of the previous section concerning 
phase stability and separation, we may now turn in more 
detail to the regime of interest for island phases. This is 
largely limited to the commensurate fillings n — 1/2, 1/3, 
and 1/4, and to the parameter range 0.1 < K/t < 0.3, 
which (Fig. 12) encompasses both the isolated phases 
which are PS endpoints, and the large-unit-cell phases. 
For n — 1/2, this region is dominated first by the flux 
phase, shown in Fig. 13. In Fig 13(a), we see the double- 
peak structure of >5'(k) with equal weigh|t-|in (0, tt) and 
(tt, O) components which is the hallmarknj of this spin 
configuration. Fig. 13(b) provides then a rare example of 
a phase where the angles between neighboring spins are 
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distributed not around the FM and AF configurations, 
but around 7r/2; our distribution is narrower than that 
in Ref. because of the larger lattice size employed. 
This spin configuration gives rise to a uniform charge 
distribution with no inhomogeneous ordering. Fig. 13(c) 
shows the real-space spin structure. 
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FIG. 13. MC results for n = 1/2 at Kit = 0.12, character- 
izing the flux phase on an 8x8 lattice, (a) Structure factor, 
(b) Angle histogram, (c) Configuration snapshot. 

By contrast, for the same filling at larger Kjt, it is 
possible to find inhomogeneous charge structures. The 
(7r/2,7r) phase of Fig. 6 exists as an endpoint both 
of phase separation and of the large-unit-cell series 
(Fig. 12). In this structure, electrons are delocalized 
across every second bond, or equivalently every FM bond 
in the 7r/2 direction, and are much more weakly present 
on the alternate AF bonds. This simple picture implies a 
stripe-like charge order with wave vector (7r,0), and the 
phase would give peaks in x-ray diffraction or electron 



microscopy experiments, which measure the charge dis- 
tribution n(r). However, because the charge density Ui is 
the same on all sites, there is no structure in the quantity 
ri(k), which is defined in Eq. (||) and readily calculated 
on a finite cluster. This situation arises only for period- 
icities 7r/2 in the spin structure factor; for all higher m 
values, n(k) and A^(k) computed from the site charges 
are indeed suitable indicators of charge order. We note 
briefiy here that by translational invariance one may in 
fact expect to find a linear superposition of equivalent 
island phases, with a uniform mean value of n^, and a 
charge order discernible only in A^(k). In the classical 
MC simulations we have shown results only for one such 
phase, which is separated by thermal barriers from its 
degenerate counterparts. 

In Fig. 14 we show a further stable configuration, which 
we call the (7r/2, 7r)-|-(7r, tt) phase. Here (Fig. 14(a)) there 
are two peaks in the structure factor, but these do not 
indicate a mixture of phases. While the histogram infor- 
mation (Fig. 14(b)) can be used only to rule out inter- 
mediate angles, it is the instantaneous MC spin configu- 
ration (Fig. 14(c)) which reveals the true nature of this 
homogeneous phase. Once again one expects a Id charge 
order for the same reasons as above. It is this phase, 
whose energy falls below the function given in Eq. (|) for 
n = 1/2 and K/t > 0.28, which breaks the large-unit- 
cell sequence, and is responsible for the wide region of 
high-contrast PS in the phase diagram of Fig. 12. 

We dwell only briefiy on the case of 3/8 filling. The 
results from the previous section show a (37r/8,7r) phase 
to be a stable endpoint in the PS regime, while the 
large-unit-cell region contains a (Stt/S, tt) member. The 
properties of these configurations are readily deduced by 
comparison with the other examples presented, and both 
have charge-ordering wave vectors of (7r/4, 0). Certain 
anomalies have been observed in experiment for filling 
n = 3/8, but these appear to be restricted to 3d systems. 
Turning to n = 1/3, the most robust island phase in 
the intermediate parameter range is (tt/S, tt), illustrated 
schematically in Fig. 1(a), and for K/t = 0.15 in Fig. 15. 
At this value of K/t, Fig. 15(a) shows a rather strong 
(7r/3, tt) peak, while Fig. 15(b) suggests a 1:2 ratio be- 
tween FM and AF bond angles despite the weak presence 
of a (7r/3, 27r/3) component. Fig. 15(c) shows the actual 
spin structure, which gives rise to a charge order at the 
wave vector (27r/3, 0), due to the higher population of 
every third site in the 7r/3 direction. This ordering is 
present in the site charge distribution function n^, which 
is shown in Fig. 16. From Fig. 16(a) it is clear that 
the charge contrast between the centre and edge sites of 
each island approaches the classical ratio of 2: 1.0 Very 
similar results are obtained for the "((1 — n)7r, tt)" phase 
(27r/3, tt) as K/t is raised beyond 0.2, as already shown 
in Fig. 7. In this state the charge-ordering wave vector 
remains (27r/3, 0). We have not been able to find a novel 
fiux phase for 1/3 filling which might be a ground state 
anywhere in the intermediate K/t regime. 
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FIG. 14. MC results for n = 1/2 at Jf/t = 0.32, character- 
izing the (7r/2, 7r)-|-(7r, tt) phase on an 8x8 lattice, (a) Struc- 
ture factor, (b) Angle histogram, (c) Configuration snapshot. 

Finally, for n = 1/4 the energy diagram (Fig. 5) in 
the region of small to intermediate Kjt is dominated by 
a single, and very robust, island phase. The extraordi- 
nary stability of the (7r/2,7r/2) phase (Fig. 1(b)) at this 
filling is clear to see by diagonalizing the 4-site square 
cluster with hopping t. This exercise yields energy lev- 
els of — 2t, 0, 0, 2t, the location of the gaps demonstrating 
immediately why the phase so favors 1/4-filling, but is 
so unfavorable at n = 1/2. Figs. 17(a)-(c) require little 
commentary, and we note only that the AF to FM angle 
ratio here is 1:1. As in Fig. 14, the charge-equivalence of 
all sites results in a homogeneous n(k) (Eq. (||)), but the 
delocalization of charge within the 2x2 squares would 
give a peak at (tt, tt) in experiments measuring ri(r). 
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FIG. 15. MC phase for n = 1/3 at Kjt = 0.15, charac- 
terizing the (7r/3,7r) phase on a 12x12 system, (a) Structure 
factor, (b) Angle histogram, (c) Configuration snapshot. 

Returning to the question of phase transitions, these 
may be considered as a function of /C/t or as a function of 
filling. In the former case the results are essentially those 
of Figs. 3-5. The only robust phases preceeding those in 
Figs. 13, 15, and 17 are FM phases, and at higher Kjt 
a short cascade of further states leads to the AF con- 
figuration. As described in Sec. Ill, the phases arising 
from MC simulations require a renormalization of their 
final energies to account for system size, and when this 
is performed the crossovers are fully consistent with the 
infinite-system results. In the experimentally more rel- 
evant case of fixed Kjt and variable filling, the results 
of Sec. IV imply that, for all but the smallest values of 
Kjt, "transitions" take the form of a differential occupa- 
tion of undoped and commensurately filled states, with 
the exception of the regime 0.2 < Kjt < 0.28 where they 
are replaced by a continuous evolution in the period of a 
large-unit-cell phase. 
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FIG. 16. Charge distribution function Ui for n = 1/3 
at K/t — 0.15, illustrating charge order of (7r/3,7r) phase 
on a 12x12 system, (a) Site charge densities: site numbers 
1-12 label the first column from bottom to top (see (b), and 
Fig. 15(c)), 13-24 the second column from bottom to top, and 
so on. (b) Charge contour plot: high densities in white, low 
in grey. 

Returning to the experiments presented in the intro- 
duction, our results justify certain, rather broad conclu- 
sions. Manganite systems which are structurally layered, 
or become layered as a result of orbital ordering in the cu- 
bic system, may indeed be susceptible to the island-phase 
phenomena, with resultant charge and spin order, dis- 
cussed here. The fundamental ingredient for this is only 
the competition between K and t intrinsic to all materials 
in the class. However, we have emphasized throughout 
the crude nature of the model we consider, and close with 
a brief discussion of the possible extensions which may 
be required to reproduce more closely the physics of real 
materials. 
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FIG. 17. MC phase for n = 1/4 at K/t = 0.12, character- 
izing the (7r/2,7r/2) phase on an 8x8 lattice, (a) Structure 
factor, (b) Angle histogram, (c) Configuration snapshot. 

One of the fundamental features of manganite systems 
is the doubly degenerate nature of the eg orbital. This 
has beea included by a number of authors, and has been 
arguedcj to be essential in accounting for the CE-typ; 
(planar in 3d) charge order observed in Lai_a;Sra;Mn03 
A further important ingredient in manganite systems is 
Jahn- Teller distortion of the local structural environment 
of each Mn ion,L3 which may act to lift the Cg orbital de- 
generacy, and also to promote charge order. Both terms 
have been included in a classical MC study of the type 
performed here,E£l albeit on very small systems. Island 
phases, in the orbital or spin degrees of freedom, were 
not among the already very rich variety of phases con- 
sidered. When two Cg orbitals pare considered, on-site 
Coulomb interactions were foundEZi to lead to the forma- 
tion of an upper Hubbard band, and to cause significant 
spectral weight shifts and broadening. Yet another term 
in many models of strongly-correlated electrons is a pos- 
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sible Coulomb repulsion between nearest-neighbor sites, 
conventionally denoted as V. This contribution acts to 
suppress phase separation, and to promote a charge or- 
dering when V competes with the hopping energy scale 
t, as noted in the Id system.B In higher dimensions, suf- 
ficiently strong V may lead to anisotropic charge order if 
the hopping is anisotropic, and more generally for weak V 
one expects a moving of phase boundaries to favor homo- 
geneous states such as the stripes and islands considered 
here. Precisely this physics was found in Ref. |2j, where 
we note that the terminology "island" phase is applied to 
mean a shrinking of the size of phase-separated regimes. 
We stress that the island phases and charge order in our 
study are intrinsic to the physics of the competing double 
exchange and superexchange, and that an additional V 
term is not required for their appearance. 

Finally, one of the major restrictions of the current 
approach is the limitation to small system sizes, which 
become smaller still on addition of the further terms 
discussed in the previous paragraph, and then still to 
largely classical considerations. The method of classical 
MC with diagonalization of the one-electron problem is 
in fact not particularly sophisticated, and we highlight 
here only two rather recent contributions which have 
the potential to reveal many more features on systems 
large enough to be considered thermodynamicallyn rep- 
resentative. These axe the variational mean- fieldEa and 
hybrid Monte Carlom techniques, both introduced for 
the double-exchange problem by the same group of au- 
thors, which allow extensions in the former case to 96^ 
systems with appropriate approximations, and in the lat- 
ter to 16^ sites with rather fewer. A last important point 
is the question of corrections to the above results due to. 
the effects of quantum fluctuations. In Id, it was foundtl 
that the boundaries between phases were moved to sig- 
nificantly larger values of K/t than predicted classically. 
While the methods presented herein do little to allow an 
assessment of fluctuation effects, these should be signifi- 
cantly smaller in 2d, both directly because of the higher 
dimensionality, and because the Id results were obtained 
with a locahzed {t2g) spin S = 1/2, whereas the classical 
limit may be no less representative of the physical situa- 
tion {S = 3/2). Thus our phase diagrams can be expected 
to be qualitatively quite accurate. One may also ask if 
quantum fluctuations would act to destroy the coherence 
of the large-unit-cell phases: because these phases are 
not spiral-ordered, and already possess the AF or FM lo- 
cal spin alignment favored by fluctuations, they may be 
assumed to be robust in this respect. 



VI. SUMMARY 

In conclusion, we have considered the possibility of 
"island" phases and associated charge order in 2d sys- 
tems, using as a model the augmented FKLM with strong 
Hund coupling. Indeed we flnd that stripe-like and island 



phases are stable at intermediate values of K/t for each 
of the commensurate fillings n = 1/2, 1/3 and 1/4. This 
result includes stability against global phase separation, 
even in the absence of additional Coulomb terms. Spi- 
ral magnetic order appears near the antiferromagnetic 
regimes at low filling or at large K/t. A variety of "flux" 
phases is possible, because the electron phase factor is 
non-trivial in all dimensions d > 1, but we find only one 
to be a stable ground state and this at n = 1/2. While 
the flux phase has a homogeneous charge distribution, 
the majority of the island phases show a charge modula- 
tion. Thus even the simple form ([|) of the FKLM repro- 
duces some of the most important experimental features 
of manganite charge and spin order. The critical values 
of K/t for transitions between ordered phases, and be- 
tween ordered and separated phases, may be identified 
rather accurately from classical considerations augment- 
ing small-system studies. 
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